Elementary processes governing the evolution of road networks 
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Urbanisation is a fundamental phenomenon whose quantitative characterisation is still inadequate. 
We report here the empirical analysis of a unique data set regarding almost 200 years of evolution 
of the road network in a large area located north of Milan (Italy). We find that urbanisation is 
characterised by the homogenisation of cell shapes, and by the stability throughout time of high- 
centrality roads which constitute the backbone of the urban structure, confirming the importance of 
historical paths. We show quantitatively that the growth of the network is governed by two elemen- 
tary processes: (i) 'densification', corresponding to an increase in the local density of roads around 
existing urban centres and (ii) 'exploration', whereby new roads trigger the spatial evolution of the 
urbanisation front. The empirical identification of such simple elementary mechanisms suggests the 
existence of general, simple properties of urbanisation and opens new directions for its modelling 
and quantitative description. 

PACS numbers: 89.75.Hc,89.75.-k,89.75.Fb,89.65.Lm 



Urbanisation is a fundamental process in human his- 
tory [1] , and is increasingly affecting our environment [2] 
and our society [3H5]. The fraction of the world popu- 
lation living in urban areas has recently grown beyond 
50% and is expected to rapidly increase in the near fu- 
ture [6]. The challenges presented by such a rapidly ur- 
banising world are complex and difficult to handle. They 
range from the increasing dependence of our society on 
fossil fuels [7] , to the emergence of socio-spatial inequal- 
ities [Slinj, and to the up-rise of serious environmental 
issues [W' . Controlling urbanisation has already proven 
to be a difficult task, and it will become even more dif- 
ficult in the near future, due to the large and unprece- 
dented magnitude of urban expansion [TTl, and to the 
intrinsic complexity of this phenomenon, which gener- 
ates cities that are continuously changing over time [H] 
in a non-homogeneous fashion [13] . Trying to understand 
the elementary spatial mechanisms that govern urbani- 
sation, leaving out specific historical, geographical, so- 
cial and cultural factors, is nowadays more important 
than ever, especially because policy makers, professionals 
and researchers are actively looking for new paradigms in 
urban planning, land management and ecology 14-fl6j . 
The existence of acultural and non-demographic drivers 
of urbanisation has been investigated using different ap- 
proaches ^17 , including classical studies in regional sci- 
ences [TS] , theory of complex systems [TOH^ , urban the- 
ory [52] and remote sensing [23]. A consistent amount 
of literature in urban history and morphology indicates 
that roads are a fundamental driver in urban evolution 
and, at the same time, one of the long-lasting constituent 



elements of urban forms [2H[2S]. In the last decade, com- 
plex networks theory has provided significant contribu- 
tions to the quantitative characterisation of urban street 
patterns [27^29] . Several studies have shown that road 
networks not only play a central role in the spatial organ- 
isation of urban areas [SO] |3T] , but are important for the 
dynamical processes occurring on them [351 I32H55] and 
for the evolution of urban systems in general |361 I37j . 
However, a quantitative analysis of the historical devel- 
opment of urbanisation in metropolitan areas is still miss- 
ing, and empirical evidence of the basic mechanisms gov- 
erning urbanisation dynamics is still lacking. Our paper 
addresses this problem, and provides a study of the evo- 
lution over two centuries of the road network in a large 
urban area located north of Milan (Italy) . 

The area under study, known as Groane, covers a sur- 
face of 125 km , includes 29 urban centres within 14 mu- 
nicipalities, and has essentially evolved along two main 
radial paths, connecting Milan to Como and Milan to 
Varese, respectively. The former path was constructed 
by the Romans during the II century B.C., while the 
latter was created during the XVI century. In the last 
two centuries, the Groane area has faced a complex pro- 
cess of conurbation, changing from a polycentric region 
into a completely urbanised area. This conurbation pro- 
cess is common to many large European metropolitan 
regions [35] • Despite local differences, it is possible to 
identify four distinct phases that characterize urbaniza- 
tion in the Groane area: i) Rural phase, (1800-1918): 
fundamentally pre-industrial, fully based on agricultural 
economy, with no major transportation infrastructures 
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FIG. 1. (a) Evolution of the road network from 1833 to 2007. For each map we show in grey all the nodes and links already 
existing in the previous snapshot of the network, and in colors the new links added in the time window under consideration. 
(b) Map showing the location of the Groane area in the metropolitan region of Milan, (c) Time evolution of the total number 
of nodes A'' in the network and of the total population in the area (obtained from census data). 



present, ii) Early-urban phase (1918-1945): between 
world wars period, still significantly based on agricul- 
ture economy, witnessing the first appearance of rail net- 
work, first small-scale sparse industrial colonisation, and 
limited expansion of rural settlements around the his- 
torical centres, iii) Urban-industrial phase (1945-1990): 
remarkable sprawled residential and industrial develop- 
ment (especially textile and mechanics) along with pop- 
ulation growth and highway construction, iv) Metropoli- 
tan post-industrial phase (1990- 2012): decline of indus- 
trial activities, slower-paced urban sprawl, former poly- 
centric organisation overwhelmed by the metropolitan 
continuum caused by the merging of expanded centres, 
increased long-range mobility due to the development of 
high speed trains and large highway systems |39j . It is im- 
portant to note that this area has never been subjected to 
overall large-scale planning efforts, one reason being that 
14 different administrative bodies preside over 14 differ- 
ent municipalities. By importing historical topographical 
and photogrammetrical data (for details see the Methods 
section) into a Geographical Information System (GIS) 
environment, we reconstructed the detailed road sys- 
tem (including minor streets) at seven different points 
in time, t = 1,2,..., 7, respectively corresponding to 
years: 1833, 1914, 1933, 1955, 1980, 1994, 2007. For each 
time, we constructed the associated primal graph [55], 
i.e. the graph where the nodes represent street junc- 
tions and the links correspond to road segments. The 
geographical information sources used to construct the 
primal graphs are reported in Table llj 



In Fig. [T^ we display the road networks at different 
times, showing how the initial small separate villages 
have grown by the addition of new nodes and links, even- 
tually merging together in an homogeneous pattern of 
streets. We denote by Gt = G{Vt,Et) the graph at time 
t, where Vt and Et are respectively the set of nodes and 
the set of links at time t. The number of nodes at time 
t is then N{t) = \Vt\, while E{t) = \Et\ is the number 
of links. By definition, we have Vt = Vt-i U AVj and 
Et = Et-i U A£'t, where AVi and Ai?t are, respectively, 
the set of new nodes and the set of new links added to 
the system in the time window ]i — 1, t]. In the following 
we will study the structure of the graph Gt at different 
times t, focusing in particular on the properties of the 
new links. As we will show, a quantitative analysis of 
the temporal evolution of the system will be able to re- 
veal in particular the existence of two simple mechanisms 
that drive the evolution of the road network over time. 



RESULTS 

Characterising the growth of the road network. 

The Groane area is basically characterised by an unin- 
terrupted growth in the period under consideration and 
displays remarkable modifications of the road system. As 
shown in Fig. [T]:, in less than two centuries, the total 
number of nodes N has increased by a factor of twenty. 



Date 
1833 
1914 
1933 
1955 
1980 
1994 
2007 



Topographical Map of Lombardy-Venetia Kingdom 

Map of Italy 

Map of Italy 

Aerial Photography Survey 

Lombardy Regional Map 

Lombardy Regional Map 

Mosaic of Urban Municipalities Plans 



Owner 

Italian Military Geographic Institute 

Italian Military Geographic Institute 

Italian Military Geographic Institute 

Italian Military Geographic Institute 

Lombardy Region 

Lombardy Region 

Lombardy Region 



Format 
Raster 
Raster 
Raster 
Raster 
Raster 
Raster 
Vectorial 



TABLE L List of geographical information sources used to construct the data set. 



from the original 255 nodes present at t — 1 (year 1833) 
to more than 5000 nodes at i = 7 (year 2007). How- 
ever, the rate of growth is not constant over time: it is 
slow from 1833 to 1933, fast from 1933 to 1980, and slow 
again from 1980 to 2007. These different growth rates 
are the signature of distinct phases of the urbanisation 
process, and are strongly related to the growth of the 
population. We indeed show that the number of nodes 
A^ is a linear function of the number of people living in 
the Groane area (see Fig. [2^) or, in other words, that the 
average number of people per road intersection in the 
area remains constant over time. 

These results highlight the peculiarity of an area where 
the population remains almost uniformly distributed in 
space, and no single urban center stands out over the 
others. In order to remove the demographic compo- 
nent and focus solely on the road network evolution, we 
adopt the number of nodes N as the natural internal 
clock of the system, and we study the change of the net- 
work properties as a function of N. The number of links 
E grows almost linearly with N (Fig. ^p, top), show- 
ing that the average node degree {k){t) — 2E{t)/N{t) is 
roughly constant over time, except for a slight increase 
from (fc) ~ 2.57 to 2.8 when going from 1914 to 1980. In 
the bottom panel of Fig. [2]d we report the total network 
length Ltot{t) — X^eeB ^(^) (where e denotes a network 
edge, and l{e) its length) which increases over time as 



N^ 



(1) 



where 7 ~ 0.54. Consequently, the average length of 
links decreases as N^~^ . This result is consistent with 
the evolution of almost regular two-dimensional lattices 
with a peaked link length distribution [221 US]- Indeed, if 
we consider a uniform distribution of nodes with density 
p on a two-dimensional surface, the typical link length 
will be of order I' ^ Xj ^^ which implies that the total 
length grows as Ltot ~ E(! ~ \/iV, thus giving a value 

7 = 1/2. 

Additional information on the structure of the road 
network can be obtained by looking at the quantity |41| 
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where N^. denotes the number of nodes of degree k (no- 
tice that we do not take into account nodes having k = 2 



in the sum, since these are not usually considered as 
proper junctions). This quantity r^v measures the rel- 
ative abundance of dead ends (corresponding to Ni) or 
T-shaped intersections (A^s), so that a small value of rj^ 
indicates a dominance of fc = 4 junctions and reveals the 
presence of a large amount of grid-like patterns. Con- 
versely, the value of r^ is closer to 1 if the network 
has numerous T-shaped crossings and dead ends. The 
plot of rN versus TV (Fig. [2j;) displays a steady decrease 
from Tfq ~ 0.87 at year 1833 to rjv — 0.835 at year 
2007, consistent with an increase of the percentage of 
four-ways crossings from N4/N ~ 11% at year 1833 to 
N4/N ~ 15.5% at year 2007. Recent studies have shown 
that the abundance of T-shaped crossings seems to be 
typical of self-organised or 'organic' urban networks, like 
those of Venice or Cairo [HTj , while the grid-like layout is 
particular to cities whose shapes are the result of large- 
scale top-down planning efforts [41], like Barcelona or 
New York. Consequently, one would conclude that what 
we are observing here is the evolution from an initial self- 
organised system to a rationally planned urban network. 
However, the Groane area at i = 1833 was actually a 
fully rural, not urban-organic network, and the succes- 
sive evolution up to year 2007 has never witnessed any 
large-scale planning whatsoever. We therefore interpret 
the result of Fig. ^ as the signature of an evolution from 
a 'pre-urban' condition, with the dominance of dead-ends 
and 3— junctions typical of rural centres in the very early 
stages of growth (i.e. still constrained by the radial con- 
vergence of major roads), to an increasingly mature 'ur- 
ban' state, in which the network expands on a relatively 
unconstrained land, further away from the dense radial 
system of original centres. At that later urban stage the 
dominance of the 4— junctions grid-like pattern does not 
result from large-scale planning, but from a piecemeal 
urbanisation, made of a multitude of small-scale, plot- 
based, scarcely coordinated developments in time, i.e. a 
substantially self-organised territorial order. 



Evolution of cells: towards homogenisation. 

Road networks are planar graphs consisting of a series 
of land cells surrounded by street segments. The statis- 
tics on the area and the shape of cells can be used to dis- 
tinguish regular lattices from heterogeneous patterns. In 
particular, it has been recently observed that for Dresden 
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FIG. 2. (a) Number of nodes N versus total population (continuous line with circles) and its linear fit (red dashed line), (b) 
Total number of links E and total network length Ltot, as a function of the number of nodes A'^. The total network length 
increases as N^'^* . (c) Value of the ratio tn between the number of nodes with degree A; = 1 and fc = 3, and the total number 
of nodes. In the inset we report the percentage of nodes having degree fc = 4 as a function of A''. Notice that the relative 
abundance of four-ways crossings increases by 5% in two centuries. 



(Germany) |30j and for a simple model of road networks 
[40] . the cell area distribution P{A) is a power law 



P{A) - A- 



(3) 



with an exponent r very close to the value 1.9. This value 
of r can be explained in terms of a lattice constructed on 
a set of nodes with density fluctuations [29 . In Fig. Isk 
we report the distribution of the cell size in the Groane 
area at i = 2007, which is indeed a power law with the 
same exponent r = 1.9 ± 0.1. However, the exponent 
changes in time, as reported in the inset: it takes a value 
r ~ 1.2 at year 1833 and converges towards r ~ 1.9 as 
the network grows. Because a larger exponent indicates 
a higher homogeneity of cell areas, we are thus witness- 
ing here a process of homogenisation of the size of cells. 
Accordingly, the relative dispersion 5a = (j{A)/iJi[A) of 
cell areas (where /i(A) and (y{A) are the average and the 
standard deviation of A, respectively) decreases from 0.5 
at year 1833 to 0.26 at year 2007, indicating that the vari- 
ance of the distribution becomes smaller as N increases. 
The diversity of the cell shapes can be quantita- 
tively characterised by the so-called shape factor $ — 
A/{'kD'^ /4)^ defined as the ratio between the area A of 
a cell and the area of the circle of diameter D circum- 
scribed to the cell [3B] . The value of the shape factor is in 
general higher for regular convex polygons, and tends to 
1 when the number of sides in the polygon increases. The 
distributions P(<&) reported in Fig.pb clearly reveal the 
existence of two different regimes: before 1933 the dis- 
tributions are single-peaked and well approximated by a 
single Gaussian function with an average of about 0.5 and 
a standard deviation of 0.25. Conversely, after 1955 the 
distributions of the shape factor display two peaks and 
can be fitted by the sum of two Gaussian functions. The 
first peak coincides roughly with the one observed before 
1933, while the second peak, centred at 0.62, signals the 
appearance after 1955 of an important fraction of regular 



shapes, such as rectangles with sides of similar lengths. 
In Fig.[3fc-d we show the cell shapes at different times (be- 
fore and after 11955) which visually confirm the findings 
of Fig. [3] The decrease in the relative dispersion 5a and 
the increase in the fraction of regular shapes suggest that 
the network undergoes an evolution towards homogeni- 
sation. This appears to result from a combined process 
that exhibits two clear patterns: i) the fragmentation of 
larger cells of natural land into smaller ones, then head- 
ing towards urbanisation by medium-large manufactures 
or services, and ii) the mostly residential urbanisation 
of peri-urban natural land in successive rings around the 
historical urban centres. In particular, this latter pat- 
tern explains the emergence of more regularly shaped 
cells after the Second World War. As urbanization took 
place increasingly further from historical main roads, res- 
idential blocks became less constrained by the triangular 
shape defined by those roads and, as a result, increas- 
ingly followed a regular rectangular shape. This regular 
shape is in fact the most efficient way of subdividing land 
for urbanization, reflecting the inner organization of the 
block into equally regular rectangular plots. This grid- 
like layout tends to be applied extensively when local 
constraints, like main roads converging into a village, do 
not force the development along different patterns (for 
example T-shaped junction dominated) [32]. 



Properties of new links: elementary processes of 
urbanisation. 



Road networks grow by the addition of new streets 
(links) and new junctions (nodes) . We focus here on the 
properties of new links by looking at their length and 
centrality. In Fig.|4]we show the cumulative distribution 
of the length of new links according to the time-section 
in which they appeared flrst. The inset shows that the 
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FIG. 3. (a) The size distribution of cell areas at f = 2007 can be fitted with a power-law p{A) ~ A^^, with an exponent 
r ~ 1.9. The inset shows the value of r at different times, (b) Distribution of cell shapes at different times, as quantified by 
the shape factor $. The shape factor of different polygons is reported at the bottom axis for comparison, (c) and (d) Maps 
showing the cell shapes (white lines) for the network as it is before 1955 (left panel) and as it is after 1955 (right panel). We see 
on the left panel that we have predominantly triangles and rectangles, while on the right panel we can observe a predominance 
of rectangles with sides of almost the same length. 



average length of new links steadily decreases over time, 
as expected from the general considerations of land frag- 
mentation reported above. More precisely, we consider, 
at each time i, the length value -^90% (i) such that 90% of 
new links at time t are shorter than £go% (i) (i.e. such that 
P{(- < ^90%) = 0.9). We notice that the value 1^^% de- 
creases in the period 1833-1933 from 625 meters down to 
325 meters, while no sensible variation is observed from 
1933 to 1994, even if the network keeps growing. In the 
last period, i.e. from 1994 to 2007, we observe another 
decrease of ^90% (i) from 325 to 225 meters. In addition, 
the relative dispersion of the length of the new links is 
almost constant and of order one, and the distribution 
does not vary too much after 1955. This phasing fits 



well with the historical development outlined in the in- 
troduction where during the rural and early-urban phase 
up to the second world war we observe the passage to 
an urban state that is then maintained along the core 
urbanisation age in the urban-industrial period until the 
1980s, followed by a different post-industrial, metropoli- 
tan regime. 

The nature of the growth process can be quantita- 
tively characterised by looking at the centrality of streets. 
Among the various centrality indices available for spatial 
networks, we use here the betweenness centrality (BC) 
143-45] , which is one of the measures of centrality com- 
monly adopted to quantify the importance of a node or a 
link in a graph. Given the graph Gt = G(Vt,Et) at time 




FIG. 4. Cumulative distributions of the length of links added 
at different times. In the inset we report the average length 
{£) of new links. 



t, the BC of a link e is defined as 



K^) = EE- 



,(e) 



(4) 






where <7ij is the number of shortest paths from node i to 
node j, while <Jij{e) is the number of such shortest paths 
which contain the link e. The quantity 5(e) essentially 
measures the number of times a link is used in the short- 
est paths connecting any pair of nodes in the network, 
and is thus a measure of the contribution of a link in the 
organisation of flows in the network. In order to evalu- 
ate the impact of a new link on the overall distribution 
of the betweenness centrality in the graph at time t, we 
first compute the average betweenness centrality of all 
the links of Gt as: 



KGt) 



1 

{N{t)-l){N{t)~-2) 



e<£Et 



(5) 



where 6(e) is the betweenness centrality of the edge e in 
the graph Gt- Then, for each link e* £ ^Et, i.e. for each 
newly added link in the time window ]t— 1, t] we consider 
the new graph obtained by removing the link e* from Gt 
and we denote this graph as Gt\{e*}. We compute again 
the average edge betweenness centrality, this time for the 
graph Gt \ {e*}. Finally, the impact Sb{e*) of edge e* on 
the betweenness centrality of the network at time t is 
defined as 



Sbie*) 



[b{Gt)-b{Gt\{e*})] 
KGt) 



(6) 



The BC impact is thus the relative variation of the graph 
average betweenness due to the removal of the link e*. 
We have measured this quantity for the different times- 
section and we report the results in Fig. [5] Remark- 
ably, the distribution of 5(,(e) displays two well-separated 



peaks (Fig. pel) . The importance of the first peak tends 
to increase in time while the second peak decreases, until 
they mostly merge into one single peak in the last time- 
section (1994-2007). In order to understand the nature 
and the evolution of the two peaks, we focus on the geo- 
graphical location of new links with different BC impact. 

In Fig. [Sk-c we report, respectively, the map of the 
network at year 1914, year 1994 and year 2007. We have 
coloured in green the links whose centrality impact be- 
longs to the left peak, and in red the links whose 5b cor- 
responds to the right peak. The new links are divided 
into two classes: green links (small (Sb, left peak) tend to 
bridge already existing streets, while red links (large ^t,, 
right peak) are usually dead-ends edges branching out of 
existing links, generating a new crossing and splitting the 
original link into two road segments. In order to better 
characterise the two classes of edges we compute fcmm(e), 
i.e. the minimum of the degrees of both endpoints of a 
new link e, and we then consider, for each value of 5b, the 
relative abundance of new links for which kmm is equal 
to 1,2,3,4, respectively. We plot the results in a gray- 
scale stacked chart (see Fig. Isfe) for the network at year 
1955 (the plots for the other networks are similar). The 
analysis of k^m confirms that the links belonging to dif- 
ferent peaks have distinct qualitative features: links in 
the right peak have kmin = 1, while links in the left peak 
have k„iin > 2. The distribution of BC impact thus sug- 
gests that the evolution of the road network is essentially 
characterised by two distinct, concurrent processes: one 
of densification (green links, left peak, i.e. lower impact 
on centrality, kmin > 2) which is responsible for the in- 
crease of the local density of the urban texture, and one 
of exploration (red links, right peak, higher impact on 
centrality, kmin = 1) which corresponds to the expan- 
sion of the network towards previously non-urbanised ar- 
eas. At a closer sight, these two patterns tend in many 
cases to appear in a temporal sequence, the former act- 
ing in preparation for the latter, i.e. exploration being 
the first phase of a cycle of urbanisation then completed 
by a second phase of densification. Obviously, since the 
amount of available land decreases over time, at earlier 
time-sections (such as in the year 1914) the fraction of 
exploration is higher, while in the 80's it becomes smaller 
until it almost disappears in 2007. The only remaining 
peak mostly corresponds to densification with new links 
having kmin > 2. 

We finally focus on another aspect of the evolution of 
street centrality, namely the relation between the age of a 
street and its centrality. In Fig. [6^ we display the links of 
the network in 2007 with a color code depending on their 
age, and in Fig. [6b we report their BC. A simple visual 
inspection shows that highly central links usually are also 
the oldest ones. In particular, the links constructed be- 
fore 1833 have a much higher centrality than those added 
at later time-sections. More precisely, the seven curves 
in panel Fig. |6j: report the cumulative distribution of the 
BC computed on the network in 2007 for the links added 
at the different times. This is defined as the probability 
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FIG. 5. The two phases of densification (grocn) and exploration (red), iUustrated for the networks at year 1914 (a), 1994 (b) 
and 2007 (c). Panel (d) shows the probability distribution of the BC impact 5b{e) for the different time snapshots. The red 
peak corresponds to exploration, and the green peak to densification. Notice that the red peak becomes smaller and smaller 
with time, and completely disappears in the last snapshot, (e) Stacked chart showing the proportion of new links with a certain 
value of kjnin. for each value of 5;, in the network at t = 1955. The peak corresponding to exploration with a high BC impact 
is entirely due to links having fc„iin ~ 1- 



P{h < x) that a link, appearing at a certain time-section, 
has a value of betweenness centrality b smaller than or 
equal to X in the final network in 2007. We notice that 
the historical structure of the oldest roads mostly coin- 
cide with the highly central links at year 2007. In par- 
ticular, the inset of Fig. [gJ: indicates that more than 90% 
of the 100 most central links in 2007 (and almost 60% of 
the top 1000) were already present in 1833. This result 
reveals that a "backbone" of highly-central routes, that 
have framed the Groane area in the rural economy of the 
pre-industrial period, has been driving the development 
of the area across two centuries of industrialisation, ur- 
banisation and de-industrialisation up to present days, 
without any major modification. 

Urban morphologists have long observed that streets 
tend to persist in time much longer than other urban ele- 
ments like buildings or land-uses |47| . Our results confirm 
that at another level, that of the topology of the street 
network, the persistence of roads over time appears to 



be due to their multifaceted impact on the urban form 
at different scales by deeply informing factors like acces- 
sibility to space and resources, land ownership and land 
values, each of these involving different social actors and 
interests. 



DISCUSSION 

In this paper we have studied the evolution, over al- 
most two centuries, of the street network in a large large 
area close to Milan (Italy). Such an area is of inter- 
est for urban studies because it displays an important 
process of urbanisation. Urbanisation is reflected in the 
growth of the road network and displays different speeds 
at different times, with a fast increase occurring between 
1933 and 1994. Our results reveal a quantitative sig- 
nature of urbanisation on the evolution of the shape and 
size of land cells, which become more homogeneously dis- 




Betweenness Centrality 

2522-10431 

1185-2521 

477-1184 

189-476 

81-188 

31-80 

0-30 



12 4 Km 

1 I I I I I I I I 




2000 4000 6000 8000 10000 



FIG. 6. Color maps indicating (a) the time of creation of each link and (b) its value of betweenness centrality (BC) at year 
2007. (c) The cumulative distribution of BC of links added at different times. The inset reports the percentage of edges added 
at a certain time which are ranked in the top n positions according to the BC. Different curves correspond to n = 100, 500, 1000. 



tributed and square-shaped. Simultaneously, we observe 
a general trend towards a larger number of 4-way junc- 
tions, as opposed to an earlier structure of dead ends 
and 3-way junctions. These structural transformations 
appear to be the result of the interplay between two con- 
current dynamics, namely densification and exploration. 
While exploration is typical of the earliest historical pe- 
riods of urbanisation, densification predominates in the 
latest. 

We were also able to quantitatively characterise the 
stability of the structure of most central streets over 
time. For instance, the most central streets in the net- 
work in 2007 largely coincide with the oldest ones. Cen- 
tral roads appear to therefore constitute a robust spatial 
backbone which remains stable over time, and charac- 
terises the evolution of the road system as a continuous 
expansion and reinforcement of pre-existing structures 
rather than as a sharp switch towards radically new con- 
figurations. The kind of evolution that we are witnessing 
in the Groane area is therefore certainly massive in over- 
all quantitative terms, reflecting unparalleled changes in 
the pace and nature of economic, technological and social 
order during the transition from a mainly agricultural to 
a post-industrial age; however, it is still an organic form 
of evolution that builds at each step on the achievements 
of the previous, confirming in the long term a resilient 
structure that tends to persist in time all over the urbani- 
sation process. These findings quantitatively support the 
hypothesis that spatial systems undergoing fine-grained 
forms of evolution tend to exhibit simple local dynam- 
ics of change, continuously expanding upon pre-existing 
structures, a morphogenetic behaviour that they appear 
to share with many living organisms and other natural 



systems |46j . 

The results of this research, including the existence of 
the two main dynamics of exploration and densification, 
together with the persistence over time of a structural 
backbone made of highly central streets, cannot be ex- 
tended to the generality of urbanisation processes unless 
supported by further investigation of cases in different ge- 
ographical and economic positions. An important point 
in our view is that this study provides quantitative results 
which can be confirmed (or falsified), a critical aspect in 
Science. 



METHODS 

Temporal Network Data. 

The area under study covers 125km^ and includes 29 
urban centers within 14 municipalities that have devel- 
oped along two main radial paths, connecting Milan to 
Como and Milan to Varese. We sampled seven street 
networks primal graphs at different times, where the 
street junctions are represented as nodes and the roads 
(or streets) are the links of the networks. We used a 
mixed procedure of ArchMap tool extensions and python 
scripts operating over a geo-database. First, the street 
network has been drawn in vectorial format on the ba- 
sis of a collection of historical areal images and historical 
maps imported in the ArcGISlO platform. Sources used 
for each temporal steps are reported in Table 1. The 
images have been imported as a .geotif extension, and 
all the layers has been projected using the coordinate 
system Monte Mario (Rome), Italy zone and Transverse 
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Mercator projection. The first geo-referencing (align- 
ment of different images) fias been performed on tlie basis 
of historical buildings and landmarks such as paths and 
roads. The redraw operation was done using an ad-hoc 
python script for creating the node layers representing 
the street junctions. The second geo-referencing opera- 
tion has been performed over the street junctions using 
the spatial analyst tool with the spatial join utilities pro- 
vided in ArchMaplO. Subsequently, for all time layers, 
using an ad-hoc python script we produced the weighted 



adjacency lists that has been used for the network analy- 
sis, the weight correspond to the real length of the street. 
We denote by Gt = G(V(, Et) the obtained primal graph 
at time i, where Vt and Et are respectively the set of 
nodes and links at time t. The number of nodes at time t 
is then N{t) — \Vt\ and the number of links is E{t) ~ \Et\. 
Using common definitions, we thus have Vt = Vt-i U AVj 
and Et = Et-iU AEt, where AVt and AEt are respec- 
tively the new street junctions and the new streets added 
in time ]t — 1, t] to the network existing at time t — 1. 
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